Spin-dependent inter- and intra-valley electron-phonon scattering in germanium 
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We investigate tlie spin-dependent electron-piionon scatterings of the L and F valleys and the 
band structure near the conduction band minima in germanium. We first construct a 16 x 16 
k ■ p Hamiltonian in the vicinity of the L point in germanium, which ensures the correctness of 
the band structure of the lowest three conduction bands and highest two valence bands. This 
Hamiltonian facilitates the analysis of the spin-related properties of the conduction electrons. We 
then demonstrate the phonon-induced electron scatterings of the L and F valleys, i.e., the intra-F/Z/ 
valley, inter-F-L valley and inter-L-L valley scatterings in germanium. The selection rules and 
complete scattering matrices for these scatterings are calculated, where the scattering matrices for 
the intra-F valley scattering, inter-F-L valley scattering and the optical-phonon and the separated 
transverse-acoustic- and longitudinal-acoustic-phonon contributions to the intra-F valley scattering 
have not been reported in the literature. The coefficients in these scattering matrices are obtained via 
the pseudo-potential calculation, which also verifies our selection rules and wave- vector dependence. 
We further discuss the Elliott- Yafet mechanisms in these electron-phonon scatterings with the k p 
eigenstates at the L and F valleys. Our investigation of these electron-phonon scatterings are 
essential for the study of the optical orientation of spin and hot-electron relaxation in germanium. 

PACS numbers: 61.72.uf, 71.70.Ej, 72.10.Di, 78.60.Fi 



I. INTRODUCTION 

The group IV materials are attractive and promising 
candidates for the achievement of spintronic devices. 
The Dyakonov-Pcrel mechanism is absent in these ma- 
terials due to the centrosymmetry and the hyper- 
fine interaction can be suppressed by isotopic purifi- 
cation, which ensures a relative long spin-dccoherence 
time.^'^'^"'^^'^"'^"^'^ Also, the silicon-based microfabrica- 
tion technology is well-developed and extensively used.^ 
Germanium (Ge), as a group IV element adjacent to sil- 
icon, shares the good spin-decoherent property and is 
fully compatible with the existing mature nanoelectronic 
technology in silicon (Si).^'^"'^^'^^ Particularly, in con- 
trast to Si, Ge shows obvious electro-optic effect as its 
direct gap (at the F point) is close to the indirect gap 
(at the L point) and lies in the infrared range. ^^"^^'^^'^^ 
Thus the optical orientation of carriers, which is free from 
the interfacial effect and the external electric and mag- 
netic fields, can be carried out effectively in Ge-based 
devices. ^^'^^'^^'^^ Moreover, compared with Si, the longer 
spin-diffusion length stemming from the larger carrier 
mobility is helpful to the spin injection and the relative 
strong spin-orbit coupling benefits the manipulation of 

spin 10,12,14,16,17 

In recent years, Ge attracts a renewed interest both ex- 
perimentally and theoretically. In the experimental side, 
Loren et al.^^-^^ and Pezzoli et al.^^ demonstrated the 
optical injection and detection of polarized electrons and 
holes in bulk Ge and Ge-based quantum wells, where elec- 
trons are pumped optically in the F valley and quickly 



scattered to the indirect valleys. In the theoretical side, 
a progress was made in investigations of the conduc- 
tion band structure and spin-dependent electron-phonon 
scattering. A compact 10x10 k ■ p Hamilto- 
nian was constructed around the L point via the method 
of invariant. ^^'^^"^^ Moreover, Tang et al.^^ derived the 
selection rules for intra- and inter-L valley electron- 
phonon scattering, and calculated the average absolute 
values of corresponding scattering elements within the 
tight-binding model. Later Li et al.^^ demonstrated the 
scattering matrices for the inter-i valley scattering and 
acoustic (AC) contribution to the intra-L valley scatter- 
ing by using the pseudo-potential method, where the ap- 
proximated wave-vector dependence of the intra-L val- 
ley scattering is derived with the combination of k p, 
pseudo-potential and group theories. Very recently Li 
et al.^^ also reported the selection rules of the intra-F- 
L valley electron-phonon scattering in the calculation of 
phonon-assisted optical transitions in Ge. However, to 
our knowledge, the complete scattering matrices of sev- 
eral important channels of the electron-phonon scattering 
have not been discussed yet, such as the intra-F valley 
scattering, the optical-phonon (OP) contribution to the 
intra-L scattering as well as the inter-F-L valley scat- 
tering. As shown in previous works, these scatterings 
arc fundamental to understand the spin dynamics in the 
optical orientation of electron spin in Ge.^°"^^'^^'^^ 

In this work, we readdress the band structure of Ge 
near the conduction band minima and study the phonon- 
induced electron scatterings in the L and F valleys. We 
first derive a spin-dependent 16x16 k p Hamiltonian in 
the vicinity of the L point in Ge from the band basis func- 
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tions at this point. ^'^"^ Compared with Ref. 18 where the 
effective electron masses are taken as granted, here we 
start from the free electron mass and straightforwardly 
obtain the renormalization of masses from this Hamilto- 
nian. It can fit the band structure of the lowest three 
conduction bands and highest two valence bands, and 
provides the eigenstates for quantitative demonstration 
of conduction electron spin properties. 

Till now, in addition to the well-known subgroup 
technique which gives the selection rules of the wave- 
vector-independent contribution to the electron-phonon 
interaction, ^^''^■^"'^^ two approaches for deriving the ex- 
plicit wave-vector dependence of scattering matrix have 
been brought iorwa.rd}^'^'^'^^'^'^ In one approach the ini- 
tial and final electronic states are expressed as the k ■ p 
eigenstates. The wave- vector dependence of scatter- 
ing matrix element is the product of that in the elec- 
tron/phonon states and crystal potential. ^^''^^ Here the 
selection rule for each analytical term is given by the 
group-theory analysis from the symmetries of the k p ba- 
sis functions, phonon state and crystal potential, while 
the coefficients of the electron-phonon scattering ma- 
trices are integrals involving the k p basis functions 
and crystal potential and can be calculated straightfor- 
wardly via the pseudo-potential method. The other ap- 
proach, i.e., the theory of invariants, utilizes the invari- 
ance of electron-phonon interaction to the symmetry op- 
erators in the corresponding space group. '^'^''^^ The in- 
variant scattering matrix consists of products of wave- 
vector-dependent irreducible tensor components and the 
bare spin-dependent matrices and takes into account the 
symmetry of phonon states. The corresponding coeffi- 
cients should be obtained via calculation with numeri- 
cal techniques such as pseudo-potential or tight-binding 
methods. Obviously, the validity of the first ap- 
proach is sensitive to the choice of k p eigenstates, and 
the second one is not limited by the k-p Hamiltonian. 
Therefore we take the invariant method, and determine 
the coefficients via fitting with our pseudo-potential cal- 
culations. 

By applying the method of invariants, we construct the 
scattering matrices and investigate the intra- and inter- 
vallcy scatterings involving both the F and four L valleys. 
The matrix elements for the intra-F and inter-F-L valley 
scattering, the OP contribution and separated transverse 
acoustic (TA) and longitudinal acoustic (LA) contribu- 
tions to the intra-L valley electron-phonon scattering are 
provided for the first time. For each phonon mode, we 
demonstrate the lowest-order wave- vector dependence of 
the scattering matrix as it is much larger than the higher- 
order terms. It should be noted that the zeroth-ordcr 
contribution to the spin-flip scattering elements does ex- 
ist in the inter- valley scattering but vanishes in the intra- 
valley case. Furthermore, with our k p Hamiltonian at 
the L point and a 14 x 14 k p Hamiltonian at the F point, 
we analyze the Elliott^® and Yafet'^^ mechanisms in these 
electron-phonon scatterings.'^^ In all the cases above, our 
calculations with pseudo-potential method confirm our 



selection rules and analytical wave-vector dependence of 
scattering matrices. 

This paper is organized as follows. In Sec. II we con- 
struct the 16 X 16 k p Hamiltonian. In Sec. HI, we in- 
vestigate the mechanisms of the phonon-induced electron 
scattering, where the intra-F/L valley scattering, inter- 
F-L valley scattering and inter-L-L valley scattering are 
studied analytically via the symmetry consideration. Be- 
sides, we calculate the scattering matrices numerically 
with pseudo-potential method. We summarize in Sec. IV. 

II. THE k p HAMILTONIAN 

Around the conduction band minimum of Ge, the 
k p Hamiltonian with the spin-orbit coupling included 
can be derived from the symmetry at the four L 
points (Va)(l,l,l), (7r/a)(-l,-l,l), {TT/a){-l,l,l) 
and (7r/a)(l, — 1, 1) with a being the lattice constant. 
At each point, the symmetry of the Bloch states is de- 
scribed by the D^d double group with the six irreducible 
representations , L^, Lg , , and Lg .^'^^''^'^ We 
first consider the vicinity of the (7r/a)(l, 1, 1) point. For 
the sake of convenience, we choose the coordinate system 
X, y, z with the z direction along the symmetry axis [111], 
so that the unit vectors of this system, related with those 
of the crystallographic frame [100], [010] and [001], i.e., 
xq, yo andzo, by i = {xQ-yo)/V2, y = (xo+yo-2io)/V6 
and z = {xo + yo + zo)/^/3. Starting from the nonrel- 
ativistic Bloch states Li, L^, L'^, , L2> in the conduction 
band and the Lg, states in the valence band and includ- 
ing the spin-orbit interaction, we obtain 16 Bloch states 
given in Table IX of Appendix A, which are used in the 
construction of the k p Hamiltonian and the analysis of 
the electron-phonon scattering. 

Due to the space inversion symmetry in bulk Ge, the 
basis functions have defined parities. For this reason, 
the superscripts and '— ' are used to represent the 
even and odd parities, respectively. It is noted that all 
the bands are 2-fold degenerate due to the time inversion 
symmetry.^ The spin-dependent perturbation Hamilto- 
nian is given by^^'^-'^ 

in which k is the electron wave vector referred to the L- 
point, the first term describes the spin-orbit interaction 
at fc = 0, me is the free electron mass, V{y) is the spin- 
independent periodic potential and tt is the generalized 
momentum operator^^ 

IT = p + 5tt , 5tv = — -—^cr x Vy(r) . (2) 

4TOeC 

The total k p Hamiltonian matrix is written in the form 
of three terms 

H = / + i?o + Hkp , (3) 
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where -ffo is the Hamiltonian matrix at the L-point and 
the iJfcp is the hnear-fc contribution describing the in- 
ter band k-Tv mixing. The diagonal components of Hq 
arc introduced in the third column of Table IX of Ap- 
pendix A. The matrix Hq also has off-diagonal com- 
ponents responsible for the interband spin-orbit mixing 
which takes place only between the band states trans- 
forming according to the equivalent spinor representa- 
tions. The nonzero off-diagonal components which stem 
from (VV^ X p) • cr terms are 

{vA\Ho\cl2) = {v3\Ho\cll) = Ai , (4a) 

(cl|iJo|c4) = (c2|i/o|c3) = A2 , (4b) 

(c9|iJo|cll) = (cl0|i7o|cl2) = A3 , (4c) 
(t;4|Fo|clO) = (z;3|i/o|c9) = -(w2|iJo|c8) 

= -(wl|iJo|c7) = A4 (4d) 

together with 10 transposed matrix elements. Here A; {I 



= 1, 2, 3, 4) are real band parameters and are listed in 
Table I. 

The linear-A; matrix can be rewritten as 



Hkp — 







(5) 



where Hcc and H^c are 12x12 and 4x12 block subma- 
trices, respectively. Taking into account that the matrix 
elements between the states of coinciding parities vanish, 
these submatrices can further be presented in the form 



Hr. 



H+- 
H-+ 



H. 



[0 ff; 



(6) 



where = {H^^'^)\ For the matrix elements oi k ■ p 

one has 



Hcc-^ik-p) 



H^c{k -p) 



I P6k+ 
Pek+ 





P2k, 


V Qik- 



-P&k^ 
Pak- 


V2P6k^ 







Q2k+ -Q2k- 
-Q2k+ -Q2k- 
-P,,k, 
-P5kz 
V2P3k^ 
-ViPsk- 

-Qi/c_ Qik+ 
P2kz -Qik- -Qik+ 
-Qik+ -P2kz 
-gifc_ -P2k, 



Pik+ 
Pik+ 
V2Pik. 






P5kz 

-Q2k+ 
-Q2k- 

Psk- 

P3k+ 

-Pifc_ 
Pik^ 


-V2Pik 





-Q2k+ 
Q2k- 
Psk- 
-P3k+ 

\ 



(7) 



(8) 



with k± — kx iky. For the matrix H~c'^, we use the 
"bra" and "ket" basis functions in the order c6 . . . cl and 
cl2 . . . c7 and, for the matrix H~^'^, the "bra" basis func- 



tions are ordered from vl to v4 and "ket" ones from c6 
to cl. For the operator k ■ Sir, the matrix elements are 
as follows 



Hcc{k ■ Stt) = 



Hyc{k ■ S-k) = 



ask^ 
aek^ 
\ploi%k^ 
2V^aQk, 


a^k^ 


(-/3i 



V 



-agk^ (/32+a5)fc+ (-/32 - "5)^- 
aefc- (-/32 + ar,)k+ {-(32 + as)^- -2^/2|32k, 
2V2aek, {(32 

-V2aek+ {(32 

a^kj^ 




a2)k+ ( 



V (/3i - a2)fc_ (-/?! - a2)/c_ 



u u 

2^a-ik, -V2a3k^ 
V2a3fc_ -2V2Q!3fc, 



-2V2(3ik, 
(3i - a2)k 



(/3i + a2)k- {—(3i — a2)k 

{(3i-a2)k- {(3i - a2)k-- 
n n 













2V2(32k, \ 


a5)k+ {P2 - "5)^+ 
a5)fc_ {-(32 + az)k^ 
a^k^ a^k^ 
a3k+ ~a3k+ ) 

a\k^ —a\k- \ 
a.\k- 

V2aifc_ -2\/2aifc^ 
2\/2aifcz \/2oL\k^ j 



(9) 



(10) 



The coefficients P\ , P2 



and ai, a2 . . . are purely imag- 



inary and Qi, Q2, /?!, (32 are real and are also listed 
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in Table I. To illustrate the applied method to cal- 
culate the matrices (7)-(10), we consider the matrix 
elements between the L-point Bloch functions '^ji = 

-^4 ; -^5 1-^6(1)' -^6(2) ^^'^ *J = ' ' ^t{l) ' ^t{2) ' 

spectively for = 1,2,3,4, formed from the single- 
group states Ly and L3 (See Table IX in Appendix A). 
Choosing the particular basis states \Ly,i') (i' = 1,2) 
and \L^,i) (i = 1,2) we calculate, by using the symme- 
try considerations, the matrix elements of operators k p 
and Ua = {h^ /2meC^) [VF(r) x k]^ between them. Let 
us denote these matrix elements by Mij and Ua-ij- Now 
if we present the states ji and $j as linear combinations 
Yji's' Gji^i-s'Cts' 1^3',*') and Y^is Cj,isas 1^3, «) , where Us 
is the spin- up state f for s = 1/2 and spin-down state 
\, for s = —1/2, then the matrix elements between the 
spinor states are given by 



i' is 



The zero point of wave vector is located at 
(7r/a)(l, 1, 1) point, a = 5.66 A is the lattice constant 
of Ge.^ One can see that the k ■ p and {WV x k) • cr 
terms couple the bases with different parities and the 
(VI/ X p) • cr terms connect bases with the same parity. 
We determine the k ■ p parameters (given in Table I) by 
fitting with an sp^dPs* tight-binding model. One can 
see from Fig. 1 that the lowest three conduction bands 
and two highest valence bands calculated from the k ■ p 
parameters fit well with the tight-binding results. 

For the other three L points, i.e., (7r/a)(— 1, — 1, 1), 
(7r/a)(— 1, 1, 1) and (7r/a)(l, — 1, 1), the basis functions 
and Hamiltonian share the same forms with those in Ap- 
pendix A and Eqs. (5)-(10), while the coordinate system 
varies. 

It is noted that a k ■ p matrix larger than or equal 
to 12 X 12 is necessary to fit the structure of the low- 
est conduction and highest valence bands without the 
remote-band influence. These bands are pertinent to the 
electron spin relaxation. Via the Lowdin partitioning,^*^ 
the effective masses of i-tli band can be expressed as 



mo 



2 Top 



E 



x{y) 



2toc l'^^J'lL,fea(fex)=o 



ft2 



(11a) 
(lib) 



TABLE I: The parameters in the 16 x 16 k ■ p Hamiltonian 
in the vicinity of the L point. These parameters are obtained 
by fitting with the band structure in an sp'^cf's* tight-binding 
model.^^ 



eV- nm 


eV 


Qi 


0.52 


Evi 


-1.119 


Q2 


0.26 




-1.365 


Pi 


0.52 z 


Eel 


0.747 


P2 


0.32 i 


Ec2 


3.990 


Pi 


0.20 i 


Ec3 


4.110 


Pi 


0.37 i 


Ec4 


8.349 


Ps 


-0.35 i 


Ecb 


8.354 


Pe 


-0.20 j 


Ecd 


9.109 


Ml 


0.03 


Ai 


0.012 




0.01 


A2 


0.089 


ai 


-O.OH 


A3 


0.041 


a2 


0.02 I 


A4 


0.097 


a-i 


0.02 i 






a4 


0.52 i 






0:5 


0.04 i 






Q6 


0.02 i 







lp_l = 



]4 1-^6(2) (^1) ^ cZ/+|.2)(-t'3)] 

j[Li i +c{L3y + tLs,) t]. 



(12a) 



(12b) 



(y) 



where L^x ^ ~zy, L^y ^ zx and ^ (X) is the spin up 
(down) eigenstate along the z direction in the correspond- 
ing coordinate system. A is the renormalization coeffi- 
cient and c is a real parameter which can be calculated 
directly by diagonalizing the k ■ p Hamiltonian. Clearly 
the conduction electron states are spin-mixing at the L 
points. In this work the Li, L^^, and L^y are taken to be 
purely real. 



For the lowest conduction band, ~ 1.36 7Tie, which 
confirms the importance of including the fourth conduc- 
tion band. 2'' 

Moreover, with our 16 x 16 fc-p Hamiltonian, it is easy 
to obtain the eigenstates of lowest conduction bands at 
the L points. 



III. ELECTRON-PHONON SCATTERING 

We investigate the electron-phonon scatterings of L 
valleys and F valley in Ge, i.e., the intra-F/L valley, 
inter-F-L valley and inter-L valley scatterings. These 
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FIG. 1: (Color online) The lowest conduction bands and high- 
est valence bands vs. wave vector near the {n/a){l, 1, 1) point. 
The origin of wave vector (the zero point in the figure) is 
taken to be at the (7r/a)(l, 1, 1) point. The green solid curves 
represent the result calculated via the k ■ p model and the 
red dashed curves describe the one obtained via the sp'^(f's* 
tight-binding (TB) model.^'^ 



scatterings are fundamental to understand the spin dy- 
namics in the optical orientation of spin and the hot- 
electron relaxation. The complete scattering 
matrices and selection rules are determined with the 
invariant method and subgroup technique. Also, we 
analyze the Elliott'^* and Yafet'^'^ mechanisms in these 
electron-phonon scatterings facilitated by the fc • p basis 
functions (See Table IX in Appendix A) and Hamilto- 
nian obtained around the L point. ■^"'"^"''■^^ The coefficients 
in the scattering matrices are obtained by the pseudo- 
potential method, which further confirms our selection 
rules and wave- vector dependence. 

We first derive the time-reversal constraint on the 
wave-vector dependence of scattering matrix. The spin- 
related scattering matrix in the centrosymmetric crystal 
can be generally written as 



(13) 



where k (k') is the wave vector of initial (final) electronic 
state and the spin eigenstates are along z direction. / is 
the 2x2 identity matrix and cr are the Pauli matrices. 
Via the time-reversal operator, it's easy to obtain the 
time-reversal constraint on the wave-vector dependence 



v4k.k' = Al,y, = A_k' 
Bk,k' = Bk k = -B_ 



,-k, 
k',-k, 



(14a) 
(14b) 



where Ak,k' (Bk,k') is purely real (imaginary). 

Moreover, the spin-orientation dependence of the 
electron-phonon scattering can be obtained easily with 
this scattering matrix. The scattering elements with 
spin eigenstates along an arbitrary direction fi = 



(sin 6 cos 4>, sin 6 sin (j), cos 9) can be expressed as 

Mk,k';tA =^k,k' + C0s6'Bk,k';z 

-I- sin6'(cos</)i3k,k';2: + sin (/)Bk,k';y), (15a) 
Mk,k';t,i = - sin6lBk,k';z + (cosO cos(f) + ism(l))Bk^k,.x 
+ (cos6'sin(/) — z cos 0)i3k,k';y- (15b) 

We further analyze the wave- vector order of scattering 
matrix element, which is given by^^''^^ 



h 



(k, s,riq|i?op|k',s' 



\ ± ^Mk,k';s, 

±1) 



(16) 



with p, V , ujq and riq representing the crystal mass den- 
sity, crystal volume, phonon frequency and phonon oc- 
cupation, respectively. The matrix element can be ex- 
pressed in order of (5q = (k' — kf) — (k — kj),^^''^^ 



M, 



k.k':s,s' 



0,s,s' ' 



-'Di^s,s'-Sq+Sci-B2,s,s'-Sq+... (17) 



with kj (kf) being the valley center of the initial (fi- 
nal) state and (5q ^ 27r/a."'^ There are vanishing (non- 
vanishing) zeroth-order contributions to the spin-flip 
scattering in the intra- (inter-) valley scatterings. Here- 
after we derive the lowest-order wave- vector dependence 
of the scattering matrix for each phonon mode only 
due to its dominant contribution compared to the corre- 
sponding higher order terms. We also give the non-zeroth 
order contributions to the spin-flip intervalley electron- 
phonon scatterings in the spherical-band-approximation 
for completeness. 



TABLE II: Phonon polarization vectors in the long- 
wavelength limit (q <^ 27r/a). Here the superscript ("— ") 
represents the in-phase (out-of-phase) . We have used the elas- 
tic continuum approximation for diamond crystal structures. 
Two TA or TO polarizations can be linearly combined into 
any other orthonormal ones. 



CtAi ' CtOi ('l) 



{qy, -qx,Q)/^/ql + 



Sa2'?to2(<i) {i^<i^,<iy<i^,-<ix -9a)/(\/gl + iJlq|) 
Ci!a.^lo(ci) 

A. Numerical method 

To obtain the coefficients of these scattering matrices, 
we derive the electron-phonon interaction and evaluate 
the matrix elements under an empirical pseudo-potential 
model by following Ref. [45] . This method has been suc- 
cessfully used in the calculation of the spin relaxation 
time,^^ the degree of circular polarization of the lumi- 
nescence across the indirect band^^ and indirect optical 
injections^^ in bulk Si, and obtains good agreements with 
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experiments. In this model, the real single particle po- 
tential V{r) = v{r — Ria) is replaced by a smooth 
pseudo potential v{r) = VL{r) + VNL{r) + Vso{r)l ■ cr 
which includes the local potential wl, the non-local one 
vnl, and the spin-orbit coupling part v^o with I being 
the orbital momentum operator. Here the subscripts i 
and a are the indices for the primitive cells of the crys- 
tal and atoms in a primitive cell separately. This pseudo 
potential is chosen to produce the same single particle 
energy as the real potential but a much smooth wave 
function around the nuclei. Therefore, the choice of the 
pseudo potential is not unique. In our calculation, it is 
taken from Ref. [48]. The calculated electron energies 
at the band edges of conduction band for F and L val- 
leys match those in the sp^d^s* tight-binding model. 
By shifting the atom position Ria by Uia and expand- 
ing V{r) with respect to Uia, the linear term in the 
expansion is the electron-phonon interaction H^p- The 
atom displacement is related to phonon operators by 
= EqA(VpWq)^^^(aqA + al^x)^'^\(^"'' ' "^hcrc q/X 
represents the phonon wave vector/mode. We calculate 
the phonon polarization vectors eqx by an adiabatic bond 
charge model. "'^ The calculated energies for the phonons 
involved in the scattering channels studied here are 10.2, 
28.6 and 33.3 meV for X3, Xi and X4 phonons, and 
7.4, 25.6, 29.3 and 35.80 meV for ig, L2', Li and Ly 
phonons, respectively. 

With the pseudo-potential method, we obtain the 
electron-phonon matrix elements for the intra-F/L val- 
ley, inter-F-L valley and inter-L valley scatterings in Ge, 
and confirm the selection rules and scattering matrices in 
all cases. The coefficients in these scattering matrices are 
determined by fitting with the pseudo-potential results 
and are listed in Tables III- VII. One can see that gener- 
ally the scattering matrix element for the spin-conserving 
process scattering is much larger (more than 50 times) 
than that for the spin-flip one. 



B. Intravalley scattering 

We first address the intra-F/i valley electron-phonon 
scatterings in Ge with the theory of invariants, where the 
LA, TA and OP contributions are aU included. ^"'^^ The 
contributions of two TA (three OP) phonon modes are 
summed up as their sound velocities (phonon energies) 
near the F point are close to each other. '^^ It should be 
noted that, beyond the symmetry of each phonon branch, 
the wave-vector dependence of polarization vectors for 
AC (belong to Fj^g irreducible representation) and OP 
(belong to FJ5 irreducible representation) phonon modes 
at the F point are well-known and compact (See Table II), 
which facilitate our derivation of the wave- vector depen- 
dence of the matrix elements. In Ge, the wave-vector- 
order analysis of the spin-flip process has already been 
performed based on the time-reversal and space- inversion 
symmetry, where the leading-order terms for the AC 
phonons are third-order [Kiq^aqn) and those for the OPs 



are second-order (Kiqm) {l,m,n £ {x, y, z}). 



1. Intra-V valley scattering 

In intra-F scattering case, we apply the conventional 
coordinate system for simplify. The Mk.k' should be an 
invariant in Oh point group. ^'^^ First, one can see that 
the zerotli-order contribution of scattering matrix van- 
ishes as 

F7®F7=F++F+. (18) 



TABLE III: The nonvanishing coefficients in the intra-F valley 
electron-phonon scattering matrix. The first two columns are 
the coefficients for the AC phonons, and the last two columns 
show the parameters in the OP case. Hi = 8.42 eV is the co- 
efficient for the LA contribution to the spin-conserving scat- 
tering. 





AC, eV-nm^ 


OP, eV- rim 




-0.61 




R2 


0.55 


4.35 


Rs 


0.60 




Ri 


-0.83 




Rs 


-1.19 


0.0055 


Re 


2.00 


0.35 



For the Oh point group, the symmetries of spin- 
dependent matrices arc given by 

/~F+ ^^F^5- (19) 

Therefore, from the invariance of the scattering matrix 
Mk,ki , the symmetries of j4k,k' and Bk.k' can be deter- 
mined by 

AkM' - rt, Bk.k' ~ r^5- (20) 

Then we can construct their explicit forms, which are 
functions of the wave vectors and phonon polarization 
vectors. It should be noted that the time-reversal con- 
straint [Eq. (14)] must be fulfilled. The symmetry- 
allowed terms are given by 

^k,k' = ^1{UXX + Uyy + U^^) (21) 

6 

Bk,k' = *E^^-Sk?k" (22) 
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with 








For D^d point group, the symmetries of the 


sp 


= (K X q)^(ua;a; + Uyj^ - 




(23a) 


dependent matrices are reduced to 


9(2) 


= (K X q)^U2^ + (K X 




(23b) 




sp 


= (K X q)ziizz, 




(23c) 




C(4) 
0^ 






(23d) 


az^L2, {cTxyCTy) ^ {Lsx, Lsy). 






(-f^'z^a: + Kxqz)Uyz, 


(23c) 




q(6) 






(23f) 





Here we omit the subscripts k, k' for simphcity in 
Eqs. (23a)-(23f). The components of Sk,k' are connected 
by the coordinate permutation. K = k' + k and q = 
k' — k = (5q. K and q are in the same order. For AC 
modes, Uafs = UlaCcj + 1P^tc,a) with a (/3) e {x,y,z}. 
In OP modes, we set Uaa = and replace UyzjUzxjUxy 
by optical vibration amplitudes Cop x^^op y^^op z- Here 
^+ and are the phonon polarization vectors in the 
long- wavelength limit, which are given in Table II. The 
coefficients and Ri are listed in Table III. 

Finally, from the numerical pseudo-potential calcula- 
tion, we can see that for the spin-conserving scattering, 
the scattering element for LA phonon branch is more 
than 2 orders of magnitude larger than that for other 
modes, as the lowest order spin-conserving element only 
exists in LA case [See Eq. (21)]. While for the spin-flip 
scattering, the OP contribution is about 4 orders of mag- 
nitude larger than the AC contribution. 



To ensure the invariance of A^k.k', the symmetries of the 
^k,k' and Bk,k' are given by 



^k,k' ~ Li, Sk,k';2 ^ L2, 

(-Bk,k';x) -Bk,k';y) ^ {L^x,L^y). 



(26) 
(27) 



Then, the explicit wave-vector of the scattering matrix 
can be constructed as well, which is much more complex 
than that in intra-F case as the symmetry is reduced. 



2. Intra-L valley scattering 

Recently, Tang et al}^ and Li et al.^^ investigated the 
intra-L electron-phonon scattering in Ge both analyti- 
cally and numerically. Tang et al}^ gave the average 
absolute values of the scattering elements for the AC con- 
tribution, and for the spin-conserving process of the OP 
contribution in the spherical-band-approximation. Li et 
al.^^ derived the approximated wave- vector dependence 
for the AC contribution, where the contributions of the 
three AC phonon modes were summed up. Here we 
demonstrate the complete and detailed wave-vector de- 
pendence for the intra-L valley scattering matrix via the 
theory of invariants, where the LA, TA and OP contri- 
butions are all considered. 

In the intra-L valley scattering, the scattering matrix 
-^k,k' should be an invariant in the D^d point group. ^'^^ 
We consider the (7r/a)(l, 1, 1) point and choose the [111] 
coordinate system for simplicity. The results in other 
three L points can be obtained by coordinate rotation. 
Here ki ~ kf = (7r/a)(l,l,l) and should be extracted 
from the initial and final wave-vectors. In this case, 
the zeroth-order contribution of the OP modes to spin- 
conserving scattering exists, 

L+®i+ = F++F+, (24) 

and the zeroth-order spin-flip elements are forbidden by 
time-reversal symmetry. -'^^ 



^k,k' = ^l{UxX + Uyy) +E2UZZ, (28) 
8 

Bk.k';. = ^EC.'^^l';., (29) 
18 

(Bk,k';a;,Bk,k';y) = Rj {Sl^},' -x^ ^^M' iy) ' (^0) 



with 



Si^^ = (uxx + Uyy)(K X q),, (31a) 

Si^^ ^Uzz{Kx(i)z, (31b) 

Si^^ = 2U^y{qxK^ - QyKy) 

- (u^x - Uyy){qxKy + qyK^), (31c) 

Si'^^ = Qziiuxx - Uyy)Kx - 2UxyKy], (SM) 

= Kz[{uxx - Uyy)qx - 2uxyqy], (31c) 
Si^^ = UyziK^qy -I- Kyq^) - u^ziKxQx - Kyqy), (31f) 

= QziUxzKy - UyzK^), (31g) 

S^^^ = Kziu^zQy - Uyzqx), (31h) 
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and 





= (Waja: +Uyy)qz{-Ky,Kx), 


(32a) 


('9(2) c(2)N 




(32b) 


Wx J ^y ) 


= (M:Ea: + Uyy)Kz{-qy, Qx), 


(32c) 






(32d) 


(^(5) C(5)) 


y^jjjj ' y y / J, ^ y y j 








(32e) 


/q{6) q{6)\ 


= Uzz{(lxKx — QyKy, —QyKx - 


qxi^y), 






(32f) 


(iS'l ^ , s^i ^ ) 

\ J.. ' t/ / 


~ Qz^ziUxx ^ ^yyi ' ^Wa^y), 


(32g) 


fc(8) c(8)\ 


— i^Qx^x Qy^y^i^xx ^yyi ~~ 








(32h) 


(5(^\5f ) 


= qzKz{-Uyz,Uxz), 


(32i) 


(5(1°) 5(1°)) 




f32i) 


(5(") 5(")) 


L\ ^a;a? ^yyj\Hx-^^x Hy y J 






— 2?7 (a K 4~ a K ) 

'^^xyydx^^y \ Hy^^xJ^ 






yujxx ^yy )\Hx^^y i Hy^^x) 






+ 2n la K — K W 

1 ^^xyy-Lx^^x Hy-'-^y)\i 


(32k) 


(5(12) 5(12)) 


y-^xz \Hx -^^y ' Hy X ) 






— 11 ^.-.i 1^ — n„ ^ 1 
^yz\Hx^^x y / "> 






- UxziqxKx - QyKy) 






- Uyz{qxKy + QyKx)], 


(321) 


(5p),5f)) 


— Kz{{'^xx ^yy^Qy '^^xyQx: 






{,^xx ^yy^Qx ^~ '^'^xyQy] J 


(32m) 


/c,(14) c-(14)x 


— Qzii'^xx ^yy)^y '^^xy-^x 






(^^xx ^yy^-^x '^^xy-^y\^ 


(32n) 


^c(15) 0(15)^ 


— ^z{^yzQy ^xzQx:^yzQx '^xzQy)^ 






(32o) 


\^x I '-'y 7 


— Qz{'^yzKy ^xz^^XT^yz^x 






(32p) 




— (K X q^zi.'i.Uxy, Uxx '^yy)^ 


(32q) 


(5(^«',Sf') 


= (K X q)ziuxz,Uyz)- 


(32r) 



Here we have again omitted the subscripts k, k' in 
Eqs. (31a)-(32r). The components of Sk,k' cannot be 
connected by the coordinate permutation now due to 
the reduced symmetry. Here Uxx + Uyy corresponds to 
Cop z fo'' OP modes, and other expressions are the same 
as those in the intra-F case. It should be noted that 
the quadratic sum of OP contributions of Ak.k' is con- 
stant [See Eq. (28)] and Table II, which is in agreement 
with the selection rule [Eq. (24)]. Moreover, Eqs. (31a), 
(31b) and (32c)-(32f) are in agreement with the approxi- 
mated analytical forms of the intra-L scattering matrix in 
Li's work.i® Obviously our wave- vector-dependent scat- 
tering matrices arc more detailed than those in previous 
works. 15'!^ 



Finally the coefficients S^, Ci and Ri in the intra-L val- 
ley electron-phonon scattering matrix are calculated from 
the pseudo-potential calculation with the non-vanishing 
ones listed in Table IV. It is interesting to see that for 
the spin-conserving scattering, unlike the intra-F case, 
the LA contribution is close to the TA one, as the lowest- 
order spin-conserving elements exist in both cases. For 
the spin-flip scattering, the ratio of OP contribution and 
AC contribution is close to that in intra-F case. Our AC 
contribution and the spin-conserving OP contribution are 
in the same order as those in previous work.^^'^^ 



TABLE IV: The nonvanishing coefficients in the intra-L val- 
ley electron-phonon scattering matrix. The first two columns 
are the coeflicients for AC phonon modes and the last two 
columns show that in OP case. 



AC 


cV 






OP 


eV/nm 


—1 


-8.6 






Hi 


43.8 


H2 


5.8 










AC 


eV-nm^ 




OP 


eV-nm 


-fti 


-3.5 X 


10" 


3 


R, 


D.U X lU 


iX2 


2.1 X 


10" 


-3 


Ro 


9 n V 1 


Ro 


4.2 X 


10" 


-3 


Rt 


1 V 1 0'"-^ 
1.0 X lU 


1X4 


0.5 X 


10" 


-4 


Rn 


D.y X lu 


Rc 


2.0 X 


10" 


-3 


Ri n 


Zr.O X ±\J 


-Re 


3.5 X 


10" 


-3 


R12 


-6.5 X 10'^ 


R7 


-1.0 X 


10" 


5 


Rl5 


2.5 X 10"^ 


Rs 


-1.2 X 


10" 


3 


R16 


-1.8 X 10"^ 


R9 


-5.9 X 


10" 


4 


Ris 


-7.0 X 10-2 


Ria 


1.0 X 


10" 


-2 






Rii 


-6.0 X 


10" 


4 






R12 


-4.9 X 


10" 


3 






Rl3 


3.3 X 


10" 


-3 






Rl4 


5.0 X 


10" 


-3 






R15 


2.0 X 


10" 


-3 






R16 


8.6 X 


10" 


-4 






Rn 


-2.5 X 


10" 


3 






Ris 


-2.0 X 


10" 


3 






Ci 


-5.6 X IQ- 


2 


C2 


1.45 X 10"^ 


C2 


-4.8 X IQ- 


2 


a 


6.0 X IQ-^ 


Cz 


1.2 X 


10" 


-1 


C7 


-8.0 X 10"^ 


C4 


6.0 X IQ- 


-3 




1.62 X 10^2 


Cs 


8.6 X 10" 


-3 






Ce 


-2.0 X 10" 


2 






C7 


-6.0 X 10" 


4 






Cg 


-2.6 X 10" 


3 
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C. Intervalley scattering 

1. Inter~T-L valley scattering 

In this part we investigate the inter-F-i valley 
electron-phonon scattering. As pointed out by Li et al.,^^ 
the zeroth-order scattering matrix elements of both the 
spin-conserving and spin- flip processes exist. Thus we 
first focus on the analysis of the zeroth-order contribu- 
tions with subgroup techniques. '^^''^^ We then investigate 
the other lowest-order non- vanishing wave- vector depen- 
dence of the scattering induced by each phonon branch 
with the method of invariants. 

Obviously the scatterings to the four L valleys are 
equivalent in unstrained bulk Ge except for the coordi- 
nate rotation. Here we consider the (7r/a)(l, 1, 1) valley 
and choose the [111] coordinate system. The lowest con- 
duction bands at the T point are basis functions of the 
irreducible representation,^ 



^7(1) — t, 



7(2: 



= ^2 



(33) 



where the spin eigcnstates are along the z direction. The 
eigenstatcs of conduction band minima at the L point 
are given in Eq. (12), which are spin-mixed. Besides, the 
phonon states at the L point, which correspond the intcr- 
T-L valley scattering, are basis functions of L^, L2', Li 
and L31 irreducible representations.^ 

The symmetry of the inter-F-L valley scattering is de- 
scribed by the D^d point group, ^'■^^ where the symmetry 
of the wave- vectors is given by 

(34) 

Here K = k'-ki-l-k-kr and q = k'-kL-(k-kr) with 
kr = (7r/a)(0,0,0) and k^ = (7r/a)(l, 1, 1). By subgroup 
techniques, we obtain a general selection rule,'^^'^'^ 



L+ ®Vj = Ly + L2' + 



(35) 



Thus only the zeroth-order contributions from the L21 
and L3' phonons exist. 

We then give a more detailed symmetry analysis on the 
zeroth-order contributions. There are various electron- 
phonon mechanisms as the L^, and Ly states are all 
two-fold degenerate. In each particular case, we perform 
the operators in Dj,d point group on the initial, final and 
phonon states. We find that many elements are forbid- 
den, and there are relations between the nonvanishing 
elements 



C2(l) 



6(2)/' 



(36) 



(rf(i)I^L3'(.,|i^(2)> -(r7"(2)l^i3'(J^6(i)>' (37) 

in which 6*2(1) stands for the rotation around the x di- 
rection for TT. ^3/(1) and £3/(2) [-^^(i) ^-nd -^^(2)] ^^'^ 
sis functions of the Ly [Lq) irreducible representation. 



which are given in Appendix A. One can see that for 
spin eigenstates along the z direction, only the L2' (^3') 
phonon branch contributes to the spin-conserving (spin- 
flip) scattering, which is in agreement with the selection 
rules in the previous work.^^ With Eq. (37) and the time- 
revesal symmetry [See Eqs. (13) and (14)], we obtain an 
additional limit on the L^i contribution 



kr.k 



r-Ki,; 

i-Bkr,ki,;y(i3' 



(38) 



Thus the inter-F-L valley scattering matrix can be ex- 
pressed as 



A 



kr.kL 



kr,k. 



XL2 



-i-Bkr,ki,;i,(i3'(2)) 
^kr,kt(i2') 



(39) 



with 

Skr,ki,;y(i3'(2)) 



-Bkr,ki;y(i3'(l)) ~ -Bkr,ki,;i/(i3')7 



(40) 

where Akr,kt(i2') [B^rMLwiLs')] is purely real (imagi- 
nary). The coefficients are given in Table V. It should 
be noted that the two non-diagonal terms come from two 
different basis functions separately, which makes the 
spin-orientation dependence of scattering matrix in this 
case different from the general form in Eq. (15). 

Now we consider the spin-orientation dependence of 
the inter-F-L valley scattering, which is given by 



kr,ki,;t,tl 



= |Mk,,k.;M(i2')l' 
= \A^rM,{L2')\^ 

4sin6'e"''^Bi 



kr,ki;t,t 



kr,ki,;i/ 



(-^3'(2)) 



--Sin0e^^i?kr,k.;,;(i3'(l)) 



I^kr,k^(i2')l^ + 7;^^""^' 



|-Bkr,ki.;y(i3')l^' 



(41) 



kr,kt;t,il 



\MkrML;1-,liB2' 



|Mk„k.;t4(i3')l^ 
2 



isin^ -e"^Bur,kLw(^3'{2)) 



kr.ki. 



;y(-^3'(i)) 



1 + cos^ e 



Bwr.MLiy{B3'W 



(42) 



where f (J,) is the spin eigenstate along an arbitrary direc- 
tion h — (sin0 cos 0, sin0 sin cos 0). One can see that 
generally both the L2' and Ly phonons are involved in 
the spin-conserving scattering and only the Ly one con- 
tributes to the spin-flip process. The spin-flip scattering 
shows obvious spin-orientation dependence, where the 
scattering is strongest (weakest) with the spin eigenstates 
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along (perpendicular to) the z direction. This anisotropy 
wiU be weakened when considering the other three chan- 
nels of inter-F-L valley scatterings, which can be ob- 
tained by the coordinate rotations. Numerically, the 
spin-conserving coefficient |^kr,ki, (^2')! is ranch larger 
(about 50 times) than the spin-flip one |i?kr,ki,;a(-^3')l 
(See Table V), and therefore the spin-conserving scatter- 
ing is nearly isotropic. 

TABLE V: The coefficients in the inter-F-L valley electron- 
phonon scattering matrix and the first- and second-order 
contributions to the spin-flip process in the spherical band- 
approximation. Akp,ki, and -Bkr.kj^;;/ represent the values 
for the zeroth-order contributions to the scattering matrix. 
Ri and Hi are the coefficients of the first-order terms in the 
scattering matrices, and Ci denote the parameters of the 
second-order ones. D\ and D2 are the parameters of the first- 
and second-order contributions to the spin-fiip process in the 
spherical-band approximation, with the spin eigenstates along 
the [111] direction. 



eV/nm 


cV 


^kr,k,,(L2') 18.21 


Ri 


0.73 


iBkr,k,^;y{L3') 0.35 


R2 


-0.34 




R3 


-2.15 X 10"^ 




Ra 


2.46 X 10"^ 




R^o 


37.64 X 10"^ 




Rg 


-38.23 X 10"^ 




R7 


1.66 




Rs 


-2.06 


cV 


cV-A 


Hi 12.0 


Ci 


12.7 


Ha -9.15 


C2 


4.54 


H3 -30.5 




6.04 


H4 20.5 




2.51 


Di(L3) 2.79 X 10"^ 


C5 


-15.5 


Di(Li) 2.14 X 10"^ 


Ce 


4.16 




C7 


7.08 






-3.98 






) 0.54 


As the zeroth-order terms 


of the 


contributions 



the -L3/L1 phonon and the spin-ffip process of the L2' 
phonon are absent, we further derive the lowest-order 
non- vanishing wave- vector dependence of these contribu- 
tions with the method of invariants. From the space- 
inversion symmetry of the electron/phonon states and 
the wave- vectors [See Eq. (34)], one finds that the L3/L1 
(L2') phonons only contribute to the odd- (even-)order 



terms in the scattering matrix. For the L3/L1 phonons, 
the first-order terms exist in both the spin-conserving 



and spin-flip processes, given by 

A{L3^) = -iiEiQy + E2Ky), (43a) 

5,(^3,) = i?iq, + i?2ifx, (43b) 
BcciLs^) = -{Rsqy + RiKy - R^q, - RqK,), (43c) 

ByiL^^) = -RsQcc - RiK,, (43d) 

A{L3y) = i(Sig^ + E2K^), (43e) 

B^iLsy) ^ RlQy + R2Ky, (43f) 

B^iLsy) - -i?3fe - RiK,, (43g) 

By{L3y) ^ RsQy + R^K y + R^Q^ + RqK^, (43h) 

and 

AiLi)^-i{E3q,+EiK,), (44a) 

B,{Li) ^ Rrqy + RsKy, (44b) 

By{Li) = -(i?7fe + RM, (44c) 



with L^x ^ —zy and L^y ^ zx being the basis functions 
of L3 irreducible representation. For the L21 phonon, the 
second-order terms are nonvanishing in both the spin- 
conserving and spin-flip processes. Nevertheless, its con- 
tribution to the spin-conserving process is negligible com- 
pared to the zeroth-order contribution Akr.kt(i2')- 
there is no zeroth-order term in the spin-flip process, the 
second-order spin-flip terms become important and read 

5,(^20 =*Ci(Kxq)„ (45) 

and 

8 

[Bx{L2^By{L2.)]=iY,C,{Si^\sl=~>), (46) 

with 



('9(2) o(2)> 
\^ X ' ^y I 


= K,{Ky,-Kx), 


(47a) 


/c(3) c(3)n 
' ^y ) 


= Kziqy,-qx), 


(47b) 


/c(4) c(4)^ 
\^x 7 ^y J 


= qz{qy,~qx), 


(47c) 


(q{5) c(5)^ 
\'-'x 7 ^y J 


= qz(Ky, -K^), 


(47d) 




^{Kl~Kl~2KxKy), 


(47e) 




= {ql - ql,-'2qxqy), 


(47f) 


\^x I '-'y J 


= {Kxqx - Kyqy, -K^qy - Kyqx)- 


(47g) 



The coefficients in Eqs. (43)-(47) arc given in Table V and 
the spin-oricntation dependence of the first- and second- 
order scattering matrices are expressed by Eq. (15). 
Note we have again omitted the subscripts k and k' in 
Eqs. (43)- (47). One can see that the first-order spin- 
conserving terms (S^) are much larger than the spin-flip 
ones {Ri). We also show the first-order contributions (of 
the L3 and Li phonons) and second-order one (of the L2' 
phonon) to the spin-flip processes of the inter-F-L scat- 
tering in the spherical-band approximation in Table V 
for completeness. 
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2. Inter-L-L valley scattering 

Recently, Tang et al}^ and Li et al}^ demonstrated 
the inter-L valley electron-phonon scattering with the 
X-point phonons (belong to X3, Xi and X4 irreducible 
representations) involved. ^'^^ Tang et al}^ gave the gen- 
eral selection rules for the zeroth-order scattering matrix 
element and calculated the average values of both the 
zeroth- and first-order contributions, where the zeroth- 
order Xs-phonon contribution is forbidden. Li et al}^ 
further derived the spin-orientation dependence and com- 
plete scattering matrices for the zeroth-order contribu- 
tions, which were expressed with coefficients Dxi.m, 
Dxi.s and Dxi,s- Here Dxi,m corresponds to the spin- 
conserving process only and the Dxi,s and Dx^.s are 
coefficients for both the spin-flip and spin-conserving 
processes. In this work we additionally derive the ex- 
plicit scattering matrix for the first-order terms of the 
Xa-phonon contribution with the method of invariants, 
which are shown to be non-negligible for the intrinsic 
electron spin relaxation in Ge in low temperature due to 
the relative low Xa-phonon energy.^* 

As the scatterings between the four L valleys are equiv- 
alent except for the coordinate rotations, we consider the 
(7r/a)(l, 1, 1) -O- (7r/a)(l, 1, — 1) case only and take the 
crystallographic frame for simplicity. The space symme- 
try of this scattering is described by one subgroup of the 
group (See Table yVj^'^A^ASM ^y^ctc the 2-fold X3 
irreducible representation in group can be deduced 
into two one-dimensional (ID) representations and the 
wave-vectors and Pauli matrices belong to ID represen- 
tations also. Based on the space symmetry of the ini- 
tial/final states, phonon states, wave-vectors and Pauli 
matrices (shown in Table VI) and the time-reversal sym- 
metry, we construct the complete scattering matrix, 

A{X^) = -tEiiq,-qy), (48a) 

5,(^3-) = R^{K, + Ky) + R^q,, (48b) 

54X3-) = i?3(9x + qy) + Ri{qx - qy) + i?5if., (48c) 

By{X^) = Rsiq, + qy) - R^iq, - qy) + R5K,. (48d) 

Here Xt ~ (sin ^ ± sin ^)(cos ^±sm^) are the 
basis functions of the X3 irreducible representation and 
the -phonon contribution is forbidden by the space 
symmetry (See Table VI). K = k' — kj + k — kf and 
q = k' - ki - (k - kf) with kj = {n/a){l, 1, 1) and kf = 
(7r/a)(l,l,— 1). The coefficients are given in Table VII 
and the spin-orientation dependence is given by Eq. (15). 
One can see that the X^ phonons contribute to both the 
spin-conserving and spin-flip processes of the inter~L-L 
scattering, and the spin-conserving one is much larger 
(about 2 orders of magnitude) than the spin-flip one. We 
also calculate the spin-flip first-order contribution of the 
X3 phonons in the spherical-band-approximation and the 
zeroth-order contributions of the Xi and X^ phonons for 
the sake of completeness (listed in Table VII) , where we 
take the same notations for the zeroth-order scattering 



matrices as those in Li's work.-'^^ The parameters are in 
the same order with those in the previous works. "'^^'^^ 



TABLE VI: The space symmetry of the final/initial elec- 
tron states, phonon states, wave-vectors and Pauli matrices 
in one subgroup of the 6*32 group for the (7r/a)(l, 1, 1) 
(7r/a)(l, 1, -1) scattering.^'^*'^'* Here a± = -j^{cr^ ± <^y), 

k± = '^i^x ± ky) and r = (a/4)(l, 1, 1). are the ba- 

sis functions of the X3 irreducible representations in the G32 
group, which are given in main text. 



g 


J 






k± 


k. 


Xt 


Li 


Lit 


(£10) 


1 


1 


1 


1 


1 


1 


1 


1 


{C2xy\T) 


1 


±1 


-1 


±1 


-1 


1 


-Lit 


Li 


iC2.y\r) 


1 


Tl 


-1 


Tl 


-1 


-1 


1 


-1 


(C2.IO) 


1 


-1 


1 


-1 


1 


-1 


Lit 


Li 


(i|r) 


1 


1 


1 


-1 


-1 


Tl 


1 


-1 


ipxy\0) 


1 


±1 


-1 


Tl 


1 


Tl 


Lit 


Li 


iPxy\Q) 


1 


Tl 


-1 


±1 


1 


±1 


1 


1 


(P.k) 


1 


-1 


1 


1 


-1 


±1 


-Lit 


Li 



TABLE VII: The coefficients in the inter-L-L valley electron- 
phonon scattering matrix. In the first two columns, the Hi 
and RiS stand for the parameters of the first-order scattering 
matrices of X3 phonons and _Di(V3) is that for the first-order 
V3-phonon spin-flip scattering element in the spherical-band- 
approximation with the spin eigenstates along the [001] di- 
rection. The last two columns show the coefficients of the 
zeroth-order contributions to the scattering matrix.^* 



cV 


c V / nm 


Hi 


6.49 




Dxi,s 


0.18 


Ri 


-3.65 X 10' 


2 




0.66 


R2 


-1.43 X 10" 


1 


Dxi.m 


6.56 


R3 


4.58 X 10" 


-2 






Ri 


-5.17 X IQ- 


2 






R5 


8.68 X 10" 


2 






Di{X3) 


6.71 X 10" 


2 
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D. Elliott and Yafet mechanisms 

In this part we analyze the Elhott'^^ and Yafet'^'^ mech- 
anisms in the spin-flip part of the electron-phonon scat- 
tering, where the electron-phonon interaction can be gen- 
erally expressed as'^^''^'^ 

h 



Uq,A-V[yo(r)/' 



(VFo(r)xp).,T]. (49) 



Here Uq is the phonon displacement vector and q (A) 
denotes the phonon wave vector (phonon mode). Vb(r)/ 
represents the bare potential and ^^^^^2 (VVb(r) x p) • cr 
stands for the spin-orbit coupling. It should be noted 
that the bare potential and the spin-orbit coupling, which 
correspond to the Elliott and Yafct processes separately, 
share the same symmetry in all the point groups of the 
crystal. Therefore their scattering matrices have the 
same analytical forms with the total electron-phonon 
scattering matrices, which are the sum of the Elliott and 
Yafet ones. 

For the spin-flip intravalley scattering, the wave- 
vector order analysis is clear in the centrosymmet- 
ric crystal. ^^''^^'^^ Our analytical results are [See 
Eqs. (23), (31) and (32)] in agreement with the previ- 
ous statements. In the AC-phonon induced scat- 
tering, the first-order {Ki) contributions to the Elliott 
and Yafet processes cancel each other completely and 
the third-order terms {Kiqmqn) remain. While in the 
OP case, the leading-order terms of both the Elliott and 
Yafet mechanisms are second-order (Kiqm) and there is 
no such perfect cancellation. Here l,m,n ^ {x,?/, z}. In 
our numerical calculation, the first-order contributions to 
the Elliott and Yafet processes in the AC spin-flip scat- 
tering are much larger (more than three orders of magni- 
tude) than the total third-order spin-flip scattering ma- 
trix elements; and both the Elliott and Yafet terms in 
OP cases are in the same order with the total spin-flip 
matrix elements. 

We then turn to the Elliott- Yafet mechanism in the 
zeroth-order spin-flip inter-P-L and inter-L valley scat- 
terings, where our pseudo-potential calculation of the 
Elliott- Yafet coefficients are listed in Table VIII. Gener- 
ally, the Elliott and Yafet contributions are in the same 
order with those of the total scattering matrices, which 
are simply sums of the corresponding Elliott- Yafet coef- 
flcients (See Tables V and VII). We further analyze the 
zeroth-order contribution of these two mechanisms with 
the P- and L-point k ■ p eigenstates [See Eqs. (33) and 
(12)], and give the exphcit expressions with the single- 
group basis functions. For the inter-P-L valley scatter- 
ing, the Elliott- Yafet matrix elements can be expressed 
as 



iBt 



,k.;,(i3') = -^(r2- t 1^4,,,, 1^3(1) t), (50) 

-c{r^t\vl,\L,^^t)l (51) 



TABLE VIII: The coefficients for the zeroth-order contribu- 
tions of the Elliott and Yafet mechanisms in the inter-F-L 
and inter-L valley electron-phonon scattering matrices. The 
superscript "E" ("Y") stands for the Elliott (Yafet) mecha- 
nism. The first two lines represent the parameters of the spin- 
flip inter-F-L scattering and the following four lines stand for 
those in the inter-L case. One can see that the coefficients in 
Tables V and VII are sums of the corresponding Elliott and 
Yafet ones in this table. 



eV/nm 



«-Bkr,ki;i;(L3') 



-0.10 
0.45 

-0.08 
0.26 
0.25 
0.41 



where f (i) is the spin eigcnstate along the [111] direc- 
tion and the superscript "E" ("Y") denotes the Elliott 
(Yafet) mechanism. One can see that the coupling be- 
tween the lowest and upper conduction bands (Li and 
L3) at the L-point is critical to the Elliott contribution, 
which is non-negligible in the total scattering matrix ele- 
ment. Here the Yafet contribution is about 4 times larger 
then the EUiott one (See Table VIII). For the inter- 
L valley scattering, we consider the scattering between 
~ a(^' ^^'^ = f-(l, 1, 1) and take the spin 

eigenstates along the [001] direction. The Elliott- Yafet 
matrix elements are given by 



:^[{L\Y \Vl\{L,y + ^coseL,,)Y) 
~{{Lly+tcoseLi,)i'\Vl\L, i%{52) 

.|J[(Llt'|V|j(L3,-fzcos0L3.)t') 
-{{Ll^+zcos9LlJV\Vl\LU')] 



mr\vi\Lu') 



+ icsin0((LL r\Vl\L, 

+ {L\Y\Vl\L,,m, (53) 
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and 



{Lir\Vl\L,,r))], (54) 
^{[^c,inei{Ll r \Vl\L,Y) 

{Lir\vi\L,, m 

-e-^^{{Lly + tcos9Ll,)i'\Vl\L, f)]}, 

(55) 



where the superscript "t" stands for the wave function at 
the point and Y (l') is the spin eigenstate along the 

[001] direction. 9 ~ arccos and (j) = j. We omit the 
terms in Eqs. (55)-(53) as c ~ 10~^. Here the mixing 
between the lowest and higher conduction bands leads to 
the Elliott process in Z?i_s(Xi) and both the Elliot and 
Yafet ones in Z?i_s(X4). The Elliott and Yafet contribu- 
tions are comparable also in this case (See Table VIII). 
The scatterings between other L points are similar to this 
case and we ignore the discussions in this work. 



IV. SUMMARY 

In summary, wc have investigated the clcctron-phonon 
scatterings of the L and F valleys, and studied the energy 
spectra near the bottom of the conduction bands with 
the L-point k ■ p Hamiltonian in Ge. We first construct a 
16x 16 k-p Hamiltonian in the vicinity of the L point with 
the double-group basis functions, of which the energy 
spectra of the lowest three conduction and highest two 
valence bands agree with the tight-binding model ones. 
The eigenstates of this Hamiltonian are useful for the 
analysis of the spin-related properties in Ge. 

We then study the phonon-induced electron scatter- 
ings of the L and F valleys in Ge, i.e., the intra-F/i val- 
ley, inter-F-L valley and inter-i valley scatterings. Via 
the symmetry consideration, we derive the selection rules 
and compact scattering matrices in all these cases, among 
which the scattering matrices for the intra-F valley, inter- 
F-L valley, OP contribution and the separated TA and 
LA contributions of the intra-L valley scatterings are ab- 
sent in the literature. We show the lowest-order wave- 
vector dependence of the scattering matrices for all the 
related phonon modes, where the zeroth-ordcr spin-flip 
scattering matrix elements are absent (present) in the 
intra- (inter-) valley cases. For completeness wc also give 
the lowest non-zcroth order contributions in the inter- 
valley cases in the spherical band approximation. The 
spin-orientation dependence of the electron-phonon scat- 
tering can be easily obtained with the corresponding scat- 
tering matrix. Our pseudo-potential calculation provides 
the coefficients of these scattering matrices, and confirms 
the selection rules and wave- vector dependence. Finally, 
we analyze the Elliott and Yafet mechanisms in these 



electron-phonon scatterings with the k ■ p eigenstates at 
the L and F valleys. This investigation provides the nec- 
essary electron-phonon scatterings for the study of the 
optical orientation of spin and the hot-electron relaxation 
in Ge. 
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Appendix A: THE k p BASIS FUNCTIONS 



TABLE IX; The basis functions representing 16 L-point Bloch 
states, 4 in the valence band and 12 in the conduction band. 



Band state 



vi) = Ljm,) 

v2)=L-{Ll,) 
vZ) 



cl) 
c2) 
c3) 
c4) 
c5) 



-'^6{i)(-^i) 
(Li) 



r + 

^6(2) 
r + 



c6)=L+(L3) 



c7) 
c8) 
c9) 



LI {LI.) 
L&(i){Ll,) 

cl0)=i6{2)(-^3' 
Cll) = ■t'6{l)(^2' 



cl2) 



M21 



{L2 



Basic functions 



^[{x - iy) i +{x + iy) 
-^{x + iy) 4.= Lg/fa) i 

Lh'id t 



1 



.X - 



- iy) t= 



-3'(1) 



St 

s l 

z{x-iy) t; i{x + iy)'' t 
—z{x + iy) I; i{x — iy)"^ 4- 
^z[(a:-iy);+(2: + ij/)t]; 
^[-(x-iyf^+(x + iyfi] 
-^z[{x-iy) i -{x + -iy) t]; 
^[{x-iy)^^+{x^iy)'' i] 
the same as L J (Lg/ ) 



see 
see 
see 

z 4, 



L^iLl,) 

Le(i)(L3') 

Le(2)iLl') 



Ho, 



E45V 

Er,. 



El 

Etc 
El, 
El, 
Etsc 



E45C 



E45C 



El 

Ee.' 
E- 



In Table IX, in order to define the transformation ma- 
trices for the 16 L-point Bloch functions used in our 
k p Hamiltonian, we explicitly give in the second col- 
umn simple examples of the corresponding basis func- 
tions. The first column presents the state notation \v,j) 
with i = 1...4 or \c,j) {j = 1...12), the double- 
group representation and, in brackets, the correspond- 
ing single-group representation. For example, the symbol 
^6(2) (-^3' ) msans the second state of the representation 
Lq originating from the valence-band representation Ly . 
The arrow f (4,) symbolizes the spin-up (down) eigen- 
state along the z direction. For each of the representa- 
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tions L^,L^ we give two different examples of the basis the basis functions, we set one of the planes ay contain 
functions. The notation of the diagonal energies Hq jj is the axes y and z in the group, 
shown in the third column. Finally, in construction of 
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